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Abstract 

In a recent paper [l| we proposed a model for calculating cross-sections of various reaction prod- 
ucts which arise from disintegration of projectile like fragment resulting from heavy ion collisions 
at intermediate or higher energy. The model has three parts: (1) abrasion, (2) disintegration of the 
hot abraded projectile like fragment (PLF) into nucleons and primary composites using a model of 
equilibrium statistical mechanics and (3) possible evaporation of hot primary composites. It was 
assumed that the PLF resulting from abrasion has one temperature T. Data suggested that while 
just one value of T seemed adequate for most cross-sections calculations, it failed when dealing 
with very peripheral collisions. We have now introduced a variable T = T{h) where h is the impact 
parameter of the collision. We argue there are data which not only show that T must be a function 
of b but, in addition, also point to an approximate value of T for a given b. We propose a very 
simple formula: T{b) = Dq + Di(As{b)/AQ) where As{b) is the mass of the abraded PLF and Aq is 
the mass of the projectile; Dq and Di are constants. Using this model we compute cross-sections 
for several collisions and compare with data. 

PACS numbers: 25.70Mn, 25.70Pq 



I. INTRODUCTION 

In a recent paper {l| we proposed a model of projectile mutifragmentation which was ap- 
plied to collisions of Ni on Be and Ta at 140 MeV/nucleon and Xe on Al at 790 MeV/nucleon. 
The model gave reasonable answers for most of the cross-sections studied. The model re- 
quires integration over impact parameter. For a given impact parameter, the part of the pro- 
jectile that does not directly overlap with the target is sheared off and defines the projectile 
like fragment (PLF). This is abrasion and appealing to the high enrgy of the beam, is calcu- 
lated using straight line geometry. The PLF has A^^ neutrons, protons and As{= Ng + Zg) 
nucleons (the corresponding quantities for the full projectile are labelled Nq, Zq and Aq). 
The abraded system A^^, Zg has a temperture. In the second stage this hot PLF expands to 
one-third of the normal nuclear density. Assuming statistical equilibrium the break up of the 
PLF at a temperature T is now calculated using the canonical thermodynamic model(CTM). 
The composites that result from this break up have the same temperature T and can evolve 
further by sequential decay (evaporation). This is computed. Cross-sections can now be 
compared with experiment. The agreements were reasonable except for very peripheral col- 
lisions and it was conjectured in l|] that the main reason for this discrepancy was due to 
the assumption of constant T over all impact parameters. 

Full details are provided in Our aim here is to improve the model by incorporating 
an impact parameter dependence of T = T{b). While we were led to this by computing the 
cross-sections of very large PLF's (which can only result from very peripheral collisions), 
the effect of temperature dependence is accentuated in other experiments. In fact these 
experiments can be used, with some aid from reasonable models, to extract "experimental" 
values for temperature T at each b. We spend considerable time studying this although our 
primary aim was and is the computation of cross-sections from a theoretical model. 



II. BASICS OF THE MODEL 



Consider the abrasion stage. The projectile hits the target. Use straight line geometry. 
We can then calculate the volume of the projectile that goes into the participant region 
(eqs.A.4.4 and A. 4. 5 of ref [2I). What remains in the PLF is V. This is a function of b. If 
the original volume of the projectile is Vq, the original number of neutrons is Nq and the 
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original number of protons is Zo, then the average of neutrons in the PLF is < Ns{b) >= 
[V{b)/Vo]No and the average number of protons is < Zs{b) >= [V{b)/Vo]ZQ] < Ns{b) > (and 
similarly < Zs{b) >) is usually a non-integer. Since in any event only an integral number of 
neutrons (and protons) can appear in a PLF we need a prescription to get integral numbers. 
Let the two nearest integers to < Ns{b) > be N^''^{b) and N^^^'ib) = N^''^{b) + 1. We 
assume that PAr^(6)=the probability that the abraded system has A^^ neutrons is zero unless 
A^,(6)is either iV™"(6) or A/'f"^(6). Let < Ns{b) >= N^^'^ib) + a where a is less than 1. 
Then P(A^f'^^(6)) = a and P{N^^''{b) = 1 - a. Similar condions apply to PzA^)- The 
probability that a PLF with A^^ neutrons and Zg protons materializes from a collision at 
impact parameter b is given by Pns,Zs{^) = Pns{^)Pzs{^)- Once this PLF is formed it will 
expand and break up into composites at a temperature T. We use CTM to obtain these. 
All the relevant details of CTM can be found in fl and 3. We will not repeat these here. 
There can be very light fragments, intermediate mass fragments (defined more precisely 
in the next section) and heavier fragments. As the fragments are at temperature T it is 
possible some of these will sequentially decay thereby changing the final population which 
is measured experimentally. Details of evaporation can be found in [l| and jj] . 



III. ARGUMENTS FOR 6-DEPENDENCE OF TEMPERATURE 

Experimental data on Mjmf as a function of Zbound (see Fig.l in jsl) probably provide the 
strongest arguments for needing an impact parameter dependence of the temperature. Here 
Mjmf is the average multiplicity of intermediate mass fragments (in this work those with 
z between 3 and 20) and ZboMnd=sum of all charges coming from PLF minus particles with 
z=l. For ease of arguments we will neglect, in this section, the difference between Zhound 
and Zs, the total charge of all particles which originate from the PLF. A large value of Zhound 
(close to Zq of the projectile) signifies that the PLF is large and the collision is peripheral 
(large b) whereas a relatively smaller value of Zbound will imply more central collision (small 
b). For equal mass collision Zbound goes from zero to Zq, the total charge of the projectile. 

The following gross features of heavy ion collisions at intermediate energy are known. If 
the excitation energy (or the temperature) of the dissociating system is low then one large 
fragment and a small number of very light fragments emerge. The average multiplicity of 
IMF is very small. As the temperature increases, very light as well as intermediate mass 
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fragments appear at the expense of the heavy fragment. The multiphcity Mjmf will grow 
as a function of temperature, will reach a peak and then begin to go down as, at a high 
temperature, only light particles are dominant. For evidence and discussion of this see |6|. 
For projectile fragmentation we are in the domain where Mjmf rises with temperature. Now 
at constant temperature, let us consider what must happen if the dissociating system grows 
bigger. We expect Mjmf will increase with the size of the dissociating system, that is, with 
Zbound- Experimental data are quite different: Mjmf initially increases, reaches a maximum 
at a particular value of Zf,ound and then goes down. 

In Fig.l we show two graphs for Mjmf, one in which the temperature is kept fixed (at 
6.73 MeV) and another in which T decreases linearly from 7.5 MeV (at 6=0)to 3 MeV at 
bmax- The calculation is qualitative. The case considered is ^^"^Sn on ^^^Sn. CTM is used 
to calculate Mjmf but evaporation is not included. Simlarly Z bound is Zg (no correction for 
z = 1 particles). Fuller calculations will be shown later but the principal effects are all in 
the graphs. Keeping the temparature fixed makes Mjmf go up all the way till Zbound=Zo is 
reached. One needs the temperature to go down to bring down the value of Mjmf as seen 
in experiment. 

IV. USE A MODEL TO EXTRACT ^^-DEPENDENCE OF TEMPERATURE 

In our model we can use an iterative technique to deduce a temperature from experi- 
mental data of Mjmf vs Zbound- Pick a 6; abrasion gives a < >. Guess a temperature 
T. A full calculation with CTM and evaporation is now done to get a Zbound and Mjmf- 
This Zbound will be close to < >. If the guessed value of temperature is too low then the 
calculated value of Mjmf will be too little for this value of Zbound when confronted with 
data. In the next iteration the temperature will be raised. If on the other hand, for the 
guess value of T, the calculated Mjmf is too high, in the next iteration the temperature 
will be lowered. Of course when we change T, calculated Zbound will also shift but this 
change is smaller and with a small number of iterations one can approximately reproduce 
an experimental pair Zbound, Mjmf- 

For the case of ^^'^Sn on ^^^Sn we provide Table I which demonstrates this. The first 
two columns are data from experiment. The next two columns are the values of Zbound and 
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Experimental 


Theoretical 


Zj-ifytiYlfi 


MjMF 


h/tq IT! H 


MjMF 


b 
(fm) 


Required T 
(MeV) 


11.0 


1.421 


11.080 


1.424 


2.912 


6.398 


15.0 


1.825 


15.094 


1.818 


3.625 


6.108 


20.0 


2.145 


19.984 


2.131 


4.4574 


5.840 


25.0 


2.010 


25.024 


2.019 


5.289 


5.520 


30.0 


1.505 


29.854 


1.545 


6.122 


5.250 


35.0 


0.920 


34.985 


0.928 


7.072 


4.970 


40.0 


0.415 


39.639 


0.424 


8.023 


4.650 


45.0 


0.193 


44.763 


0.196 


9.331 


4.350 


47.0 


0.156 


46.512 


0.154 


9.925 


4.260 


49.0 


0.135 


48.425 


0.130 


10.876 


4.190 



TABLE I: Best fit and experimental values for Sn on ^Sn. The first two columns are data from 
experiment. The next two columns are the values of Zbound and Mjmf we get from our iterative 
procedure. These values are taken to be close enough to the experimental pair. These are obtained 
for a value of h (fifth column) and a temperature T (sixth column). 

MiMF we get from our iterative procedure. These values are taken to be close enough to the 
experimental pair. These are obtained for a value of h (sixth column) and a temperature T 
(fifth column). Table II provides similar compilation for ^°^Sn on ^^^Sn. 
Having deduced once for all such "experimental data" of T vs b, one can try simple parametri- 
sation hke T(b) = Cq + Ci * b + C2 *b'^.... and see how well they fit the data. We show this 
for the two cases in Fig. 2. 

In Fig. 3 using such parametrised versions of T we compute Mimf vs Zbound and compare 
with experimental data. Except for fluctuations in the values of Mimf for very low values 
of Zbound the fits are very good. We will return to the cases of fluctuations in a later section. 
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Experimental 


Theoretical 


'Zi l)oimd 


MiMF 


Z })OUTld 




b 
(fm) 


Required T 
(MeV) 


15.0 


1.690 


14.816 


1.583 


3.886 


6.200 


20.0 


1.923 


19.865 


1.906 


4.698 


5.740 


21.0 


1.984 


21.207 


1.976 


4.930 


5.705 


25.0 


1.749 


24.913 


1.758 


5.510 


5.320 


30.0 


1.079 


30.356 


1.075 


6.438 


4.900 


35.0 


0.581 


35.252 


0.602 


7.366 


4.600 


40.0 


0.223 


40.123 


0.225 


8.410 


4.210 


45.0 


0.201 


44.676 


0.199 


9.802 


4.100 


47.0 


0.201 


47024 


0.159 


10.876 


4.000 



TABLE II: Same as Table 1, except that here the projectile is ^'^'Sn instead of Sn. 
V. TEMPERATURES EXTRACTED FROM ISOTOPE POPULATIONS 



In the preceding sections we have extracted temperatures T (combining data and model) 
at values of b (equivalently at values of Zbound)- This is a new method for extracting temper- 
ature. A more standard way of extracting temperatures is the Albergo formula which has 
been widely used in the past (for a review see, for example, 10, 11|). In [[3], Figs. 24 and 25] 



temperatures at selected values of Zbound/ Zq were extracted from populations in p''^IIe,^'''Li] 
and ['^'^Be,^'^Li] using Albergo formula. These temperatures are compared in Fig.4 with a 
typical temperature profile deduced here. It is gratifying to see that such different methods 
of extraction still give reasonable agreement. 



VI. FLUCTUATIONS IN Mjmf FOR SMALL Zbound 

For small values of Zbound the measured Mjmf shows considerable fluctuations as we go 
from one value of Zbound to another (see Fig.3). Our model does not reproduce these although 
general shapes are correct. Statistical models are not expected to show such fluctuations 
but let us get into some details which (a) give a clue how such fluctuations may arise and 
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(b) why our model misses them. The reader who is not interested in such details can skip 
the rest of this section without loss of continuity. 

For definiteness, consider the case of ^^''Sn on ^^^Sn. By the definition of IMF, {z > 2), 
MjMF is when Z hound is 2. Consider now Zhound = 3. The most direct way one can have 
this is if the PLF has Zs='i. Taking very simplistic point of view, suppose, this also has 
Ns = 3, that is, the PLF is ^Li. This is stable and we immediately get Mjmf = 1- This is 
indeed the experimental value. The case of Zbound=4: may arise if the PLF is Zs = 4, A^^ = 4, 
i.e., if the PLF is ^Be. But ^Be is unbound and will break up into two a particles which are 
not IMF's. Thus Mjmf drops to zero. In experiment this falls to about 0.3 rather than 0. 
The simple fact that ^Li (and an excited state of ^Li) is particle stable whereas the states 
of '^Be are not is not embedded in our liquid-drop model for ground state and Fermi-gas 
model for excited states. Our description gets better for larger nuclear systems but for very 
small systems quantum mechanics of nuclear forces causes rapid changes in properties as 
one goes from one excited state to another and one nucleus to another. Our model can not 
accommodate this. 

Let us go back to the case of Zs=3 and treat it more realistically. Using the abraison 
model, when Zg is 3, PLF can have A^s=3,4 and 5. Probabilities for higher and lower values 
of Ng are small. Following our model we get Mjmf ~ 0.94 with Z^ound slightly less than 3. 
When Zg is 4, significant probabilities occur forA^s=5,6 and 7. Following our model we get 
a small increase in Mjmf with Z bound whereas experimentally Mjmf falls. This discrepancy 
happens because in the Fermi-gas model there is very little difference between properties 
of ground and excited states of Li and Be whereas, in reality they are very different. A 
much more ambitious calculation for very small dissociating systems with Zg between 3 and 
7 where we take binding energies and values of excited state energies from experiments (this 
becomes more and more unwieldy as Zg increases) is under way. 

Fig. 7 in ref [3| shows that SMM calculations are able to reproduce the fluctuations faith- 
fully. Actually in those calculations the occurrences of Zg,Ng with associated are not 
calculated but guessed so that the ensemble produces the data as_ faithfully as possible. For 
further details how these calculations were done please refer to 
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VII. TOWARDS A UNIVERSAL TEMPERATURE PROFILE 



Knowing the temperature profile T = T{b) in one case, say ^^^Sn on ^^^Sn, can we 
anticipate what T = T{b) will be like in another case, say, ^^Ni on ^Be ? In both the cases 
bmin is zero and bmax is R1+R2 yet we can not expect the same functional form T = T{b/bmax) 
for both the cases. In the first case, near 6 = a small change in b causes a large fractional 
change in the mass of the PLF whereas, for ^^Ni on ^Be, near 6=0, a small change in b 
causes very little change in the mass of the PLF. Thus we might expect the temperature to 
change more rapidly in the first case near 6=0, whereas, in the second case, the temperature 
may change very little since not much changed when 6 changed a little. In fact, for Ni on 
Be, transport model calculations, HIPSE (Heavy Ion Phase Space Exploration) and AMD 
(Antisymmetrised Molecular Dynamics) find that starting from 6=0, excitation energy/per 
particle changes very little in the beginning [l^. In terms of our model, this would mean 
that for Ni on Be, T would be slow to change in the beginning. 

We might argue that a measure of the wound that the projectile suffers in a heavy ion 
collision is 1.0 — As/Aq and that the temperature depends upon the wound. Thus we should 
expect T = T{As{b)/AQ). Just as we can write T(6) = Cq + Ci * 6 + C2 * 6^ + ... so also we 
could expand in powers of As{b)/Ao,i.e., T{b) = Dq + Di{As{b)/Ao) + D2{As{b)/Aof + ... 
We try such fits to the "experimental" temperature profile given in Tables I and II. From 6 
we deduce As{b)/AQ and plot T as a function of As{b)/AQ. A linear fit appears to be good 
enough (Fig. 5). 

The specification that T(6) = Dq + Di{As{b) / Aq) has profound consequences. This means 
the temperature profile T{b/b,nax) of ^^^Sn on ^^^Sn is very different from that of ^^Ni on 
^Be. In the first case As{b)/AQ is nearly zero for 6 = 6mm=0 whereas in the latter case 
As{b)/Ao is ^ 0.6 for 6 = 6^i„=0. For Do=7.5 MeV and Di=-4.5 MeV, the temperature 
profiles are compared in Fig.6. Even more remarkable feature is that the temperature profile 
of ^^Ni on ^Be is so different from the temperature profile of ^^Ni on ^^^Ta. In the latter 
case bmin = Rxa — Rm and beyond bmm, As{b)/Ao grows from zero to 1 for bmax- This is 
very similar to the temperature profile of ^^^Sn on ^^^Sn. 
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VIII. FORMULAE FOR CROSS-SECTIONS 



Now that we have estabhshed that temperature T should be considered impact parameter 
h dependent, let us write down how cross-sections should be evaluated. We first start with 
abrasion cross-section. In eq.(l) of ^j], the abrasion cross-section was written as 

aa,N.,z. = 27T J bdhP^Mb) (1) 

where Pns,zA^) is the probability that a PLF with A^^ neutrons and Zg protons emerges 
in collision at impact parameter b. Actually there is an extra parameter that needs to be 
specified. The complete labelling is aa,Ns,Zs,T if we assume that irrespective of the value of 
6, the PLF has a temperature T. Here we have broadened this to the more general case 
where the temperature is dependent on the impact parameter h. Thus the PLF with Ng 
neutrons and Zs protons will be formed in a small range of temperature (as the production 
of a particular N^, Zg occurs in a small range of h). 

To proceed, let us discretize. We divide the interval hmin to hmax into small segments of 
length A6. Let the mid-point of the i-th bin be < 6j > and the temperature for collision at 
< 6i > be Tj. Then 

(ya,N^,Zs =^<^a,Ns,Zs,Ti (2) 
i 

where 

(^a,Ns,Zs,% = 27r < hi > AbPNs,zA< h >) (3) 

PLF's with the same Ns, Z^ but different Tj's are treated independently. The rest of the 
calculation proceeds as in l|. If, after abrasion, we have, a system Ng^ Zg at temperature T,, 
CTM allows us to compute the average population of the composite with neutron number n, 
proton number z when this system breaks up (this composite is at temperature Tj). Denote 
this by M^i'^'''^\ It then follows, summing over all the abraded Ng, Zg that can yield n, z 
the primary cross-section for n, z is 



^ n,z 



E M:^f-^^a.,^^,z.,T, (4) 



N,,Zs,Ti 

Finally, evaporation from these composites n, z at temperatures Tj is considered before 
comparing with experimental data. 
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IX. CROSS-SECTIONS FOR DIFFERENT REACTIONS 



We will now show some results for cross-sections using our model and compare with 
experimental data. We first show results for ^^^Sn on ^^^Sn and ^'^''Sn on ^^^Sn at 600 
MeV/nucleon beam energy. The experimental data are plotted in 0] and the data were 
given to us, thanks to Prof. Trautmann. The differential charge distributions and isotopic 
distributions for ^°''Sn on ^^^Sn and ^^^Sn on ^^^Sn were theoretically calculated using T{b) = 
Co + Cih and and also T{h) = Dq + Di{As{h) / Aq). So long as the temperature values at 
the two end points of h are the same, the answers did not differ much. In Fig. 7 we have 
shown results for T varying linearly with h with T^ax = 7.5MeV and Tmm=3 MeV. At each 
Zbound, the charge distribution and isotopic distributions are calculated separately and finally 
integrated over different Zbound ranges. The differential charge distributions are shown in 
Fig. 7 for different intervals of Zbound/ Zq ranging between 0.0 to 0.2, 0.2 to 0.4, 0.4 to 0.6, 
0.6 to 0.8 and 0.8 to 1.0. For the sake of clarity the distributions are normalized with 
different multiplicative factors. At peripheral collisions (i.e. 0.8<Zbound/Zo<1.0) due to 
small temperature of the projectile spectator, it breaks into one large fragment and small 
number of light fragments, hence the charge distribution shows U type nature. But with the 
decrease of impact parameter the temperature increases, the projectile spectator breaks into 
large number of fragments and the charge distributions become steeper. In Figs. 8 and 9 the 
integrated isotopic distributions over the range 0.2<Zbound/Zo<0.8 for Beryllium, Carbon, 
Oxygen and Neon are plotted and compared with the experimental result for ^°^Sn on ^^^Sn 
and ^^^Sn on ^^^Sn reaction respectively. 

Rest of the cross-sections shown all use T{b) = 7.5MeV-{As{b)/AQ)4.5Mey. First, in 
Fig. 10 the calculations of Fig. 7 are redone but with the above parametrization. Next we 
look at data for ^^Ni on ^Be and ^^^Ta at beam energy 140 MeV/nucleon done at Michigan 
State University. The data were made available to us by Dr.Mocko (Mocko, Ph.D. thesis). 
Calculations were also done with ^^Ni as beam. Those results agree with experiment equally 
well but are not shown here for brevity. The results for ^^Ni on ^Be and ^^Ni on ^^^Ta are 
shown in Figs. 11 to 14. The experimental data are from 13|]. The chief difference from 
results shown in jl| is that we are able to include data for very peripheral collisions. Next 
we look at some older data from ^^^Xe on ^^Al at 790 MeV/nucleon [l^. Results are given 
in Figs. 15 and 16. 
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The parametrization T{b) — 7.5MeV-{As{b)/Ao)4:.5MeV was arrived at by trying to fit 
many reaction cross-section data. Better fit for Mjmf vs. Zhound for Sn isotopes is found 
with shghtly different values: T{b) = 7.2MeV-(A,(6)/Ao)3.2MeV. 

X. SUMMARY AND DISCUSSION 

We have shown that there are specific experimental data in projectile fragmentation which 
clearly establish the need to introduce an impact parameter dependence of temperature T 
in the PLF formed. Combining data and a model one can establish approximate values of 
T — T{h). The model for cross-sections has been extended to incorporate this temperature 
variation. This has allowed us to investigate more peripheral collisions. In addition, the 
impact parameter dependence of temperature appears to be very simple: T{h) = Dq + 
Di{As{h) / Aq) where Dq and Di are constants, Asip) is the mass of the PLF and Aq is the 
mass of the projectile. With this model, we plan to embark upon an exhaustive study of 
available data on projectile fragmentation. 
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FIG. 2: (Color Online) Impact parameter dependence of temperature for ^''^Sn on ^^^Sn (left 
panels) and ^^^Sn on ^^^Sn reactions (right panels). The red squires in the upper panels represent 
the extracted temperatures (sixth column of table-I and II) and the blue dotted lines are linearly 
decreasing temperature profile from 7.5 MeV to 3 MeV. The blue dotted lines of middle and 
lower panels represent fitting of extracted temperatures (red squires) with T{b) = Cq + Ci * b and 
r(6) = Co + Ci * 6 + C2 * 62 equation respectively. The unit of Cq is MeV, Ci is MeVfm"! and C2 
is MeVfm~2. 
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FIG. 3: (Color Online) Mean multiplicity of intermediate-mass fragments Mjmf, as a function 
of Zbound for ^'^''Sn on ^^^Sn (left panel) and ^^^Sn on ^^^Sn (right panel) reaction calculated using 
linearly decreasing temperature from 7.5 MeV to 3 MeV (red solid lines) and T(b) = Co + Ci*b + 
C2 * 6^ profile (blue dotted lines). The experimental results are shown by the black dashed lines. 
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FIG. 4: (Color Online) Comparison of theoretically used temperature profiles (i) temperature 
decreasing linearly with impact parameter from 7.5 MeV to 3 MeV (red solid lines), (ii)T(6) = 
Co + Ci * b + C2 * fitting temperature (blue dotted lines) with that deduced by Albergo formula 
from experimental data (black points with error bars) for ^'^''Sn on ^^^Sn (left panel) and ^^^Sn on 
^^^Sn (right panel). 
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FIG. 5: (Color Online) Fitting of extracted temperatures (red squires) with T{b) = Dq + 
Di{As{b)/Ao) (blue dotted lines in upper panels) and T{b) = Dq + Di{As{b) /Aq) + D2{As{b) / Aq)'^ 
profile (blue dotted lines in lower panels) for ^^''Sn on ^^^Sn (left panels) and ^^'^Sn on ^^^Sn (right 
panels). The units of Dq, Di and D2 are MeV. 
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FIG. 6: (Color Online) Temperature profile for ^^Ni on ^Be (black solid line), ^^Ni on ^^^Ta 
dotted line) and ^^''Sn on ^^^Sn (blue dashed line) by considering T = 7.5 — 4.5(^s(6)/Ao) 
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Proton Number(Z) 

FIG. 7: (Color Online) Theoretical total charge cross-section distribution (red solid lines) for 
-^°^Sn on ^^^Sn (left panel) and ^^^Sn on ^^^Sn reaction (right panel) sorted into five intervals of 
Zbound/ZQ ranging between 0.0 to 0.2, 0.2 to 0.4, 0.4 to 0.6, 0.6 to 0.8 and 0.8 to 1.0 with different 
multiplicative factors 10~^, 10~^, lO'', 10^, 10^ respectively. The experimental data are shown by 
black dashed lines. Theoretical calculation is done using linearly decreasing temperature from 7.5 
MeV at 6=0 to 3 MeV at hmax- 
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FIG. 8: (Color Online) Theoretical isotopic cross-section distribution (circles joined by dashed 
lines) for ^'^'^Sn on ^^^Sn reaction summed over 0.2<Zhound/Zo<0.8. The experimental data are 
shown by black squires. Theoretical calculation is done using linearly decreasing temperature from 
7.5 MeV at 6=0 to 3 MeV at bmax- 
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FIG. 9: (Color Online) Same as Fig. 8, except that here the projectile is ^^"^Sn instead of ^"''Sn. 
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FIG. 10: (Color Online) Same as Fig. 7 except that here the temperature profile is T(b) = 
7.5MeV — {As(b)/Ao)4:.5MeV instead of linearly decreasing temperature from 7.5 MeV at b=0 to 

3 MeV at bmax 
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FIG. 11: (Color Online) Total mass (left panel) and total charge (right panel) cross-section 
distribution for the ^^Ni on ^Be reaction. The left panel shows the cross-sections as a function 
of the mass number, while the right panel displays the cross-sections as a function of the proton 
number. The theoretical calculation is done using temperature decreasing linearly with As/Aq 
from 7.5 MeV to 3.0 MeV (dashed line) and compared with the experimental data (solid line). 
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FIG. 12: (Color Online) (Color Online) Same as Fig. 11 except that here the target is 
instead of ^Be. 
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FIG. 13: (Color Online) Theoretical isotopic cross-section distribution (circles joined by dashed 
lines) for ^^Ni on ^Be reaction compared with experimental data (squares with error bars). 
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FIG. 14: (Color Online) Same as Fig. 13 except that here the target is ^^^Ta instead of ^Be. 
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FIG. 15: (Color Online) Total charge cross-section distribution for the ^^^Xe on ^^Al reaction. The 
theoretical calculation is done using temperature decreasing linearly with 74^/^0 from 7.5 MeV to 
3.0 MeV (dashed line) and compared with the experimental data (solid line). 
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FIG. 16: (Color Online) Theoretical isotopic cross-section distribution (circles joined by dashed 
lines) for ^^^Xe on ^''Al reaction compared with experimental data (squares with error bars). 
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